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ABSTRACT 

r ~~ ] Tidal dissipation in compact white dwarf (WD) binary systems significantly in- 

Ph flucnces the physical conditions (such as surface temperature and rotation rate) of the 

WDs prior to mass transfer or merger. In these systems, the dominant tidal effects in- 
volve the excitation of gravity waves and their dissipation in the outer envelope of the 
star. We calculate the amplitude of tidally excited gravity waves in low-mass (0.3Mq) 
q helium-core (He) WDs as a function of the tidal forcing frequency w. Like carbon- 

S-) oxygen (CO) WDs studied in our previous paper, we find that the dimensionless tidal 

torque FiiS) (inversely proportional to the effective tidal quality factor) depends on u> 
i ^ i in an erratic way. On average, F (cj) scales approximately as w 6 , and is several orders 

of magnitude smaller for He WDs than for CO WDs. We find that tidal torques can 
begin to synchronize the WD rotation when the orbital period is less than about a 
^ hour, although a nearly constant asynchronization is maintained even at small periods. 

We examine where the tidally excited gravity waves experience non-linear breaking 
or resonant absorption at a critical layer, allowing us to estimate the location and 
magnitude of tidal heating in the WD envelope. We then incorporate tidal heating in 
the MESA stellar evolution code, calculating the physical conditions of the WD as a 
y—i function of orbital period for different WD models. We find that tidal heating makes 

CN a significant contribution to the WD luminosity for short-period (~ 10 min) systems 

such as SDSS J0651+2844. We also find that for WDs containing a hydrogen envelope, 
tidal heating can trigger runaway hydrogen shell burning, leading to a nova-like event 
V~j before the onset of mass transfer. 
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1 INTRODUCTION 

In the last decade, compact white dwarf (WD) binaries (with orbital periods in the range of minutes to hours) have become 
increasingly important for several topics in astrophysics. The orbits of these systems decay via the emission of gravitational 
waves, constituting the largest signals for next generation space-based gravitational wave detectors. Systems of sufficiently 
short orbital period will merge within a Hubble time, the result of which may create a variety of exotic astrophysical systems, 
e.g., isolated sdB and sdO stars, R CrB stars, AM CVn binaries, or high-mass neutron stars. Of most importance, merging WDs 
may trigger type la supernovae (e.g., Webbink 1984; Iben & Tutukov 1984). Recent observations and numerical simulations 
have provided some support for such "double degenerate" progenitors of SNe la. (e.g., Gilfanov & Bogdan 2010; Di Stefano 
2010; Maoz et al. 2010; Li et al. 2011; Bloom et al. 2012; Schaefer & Pagnotta 2012, Gonzalez Hernandez et al. 2012). 

The outcome of a WD binary merger depends on the masses and compositions of the WDs and their pre-merger conditions 
(e.g., Segretain et al. 1997; Yoon et al. 2007; Loren-Aguilar et al. 2009; van Kerkwijk et al. 2010; Dan et al. 2012; Raskin et 
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al. 2012). Despite the broad significance of WD mergers in astrophysics, detailed studies of the pre-merger conditions have 
been relatively scarce. Most studies have focused on equilibrium (non-dynamical) tides (e.g., Iben et al. 1998, Willems et 
al. 2010) or have parameterized the tidal effects (e.g., Piro 2011). None of these studies have sought to calculate both the 
magnitude and location of tidal heating, and none of them can be used to predict observational signatures of tidal heating. 
Such predictions are becomingly increasingly important as ongoing surveys continue to uncover new compact WD binary 
systems (e.g. Mullally et al. 2009; Kulkarni & van Kerkwijk 2010; Steinfadt et al. 2010a; Kilic et al. 2012; Brown et al. 2011; 
see Marsh 2011 for a review). 

This paper is the fourth in a series (see Fuller & Lai 2011, 2012a, 2012b, hereafter Papers I, II, and III) where we system- 
atically study dynamical tides and their observational consequences in compact WD binaries. In Paper I, we calculated the 
tidal excitation of discrete g-modes in carbon-oxygen (CO) WDs. We showed that the excited g-modes reach very non-linear 
amplitudes near the surface of the star, even far from resonance. Thus, rather than exciting discrete g-modes, the binary 
companion will excite a continuous train of gravity waves that propagate towards the surface of the WD, where they are likely 
dissipated through a combination of non-linear processes and radiative damping. 

In Paper II, we calculated the gravity wave amplitude as a function of orbital frequency using an outgoing wave boundary 
condition (i.e., assuming the waves damp completely in the outer layers of the WD rather than reflecting at the surface). Our 
calculations showed that the train of gravity waves is launched at the C-He composition gradient in the CO WD models we 
used. We then computed the energy and angular momentum flux carried by the waves in order to predict the orbital and 
spin evolution of WDs in compact systems. We found that tidal effects are negligible at large separations (orbital periods 
larger than about an hour), but become increasingly important at smaller periods, causing the WDs to be nearly synchronized 
upon merger. Furthermore, we found that the heating rate can be quite large at short orbital periods (exceeding 100L Q just 
before merger, depending on the system parameters), potentially leading to an observable tidal heating signature. In Paper 
III, we showed that localized tidal heating in the semi-degenerate region of the hydrogen envelope can lead to thermonuclear 
runaways, "tidal novae", potentially burning the hydrogen layer off the WD in an explosive event similar to a classical nova. 

In this paper, we extend our calculations to models of low mass (M « O.SMq) helium (He) WDs, which are relatively 
common among observed short-period WD binary systems (Kilic et al. 2012). Since these WDs do not have a C-He composi- 
tional gradient (although they do contain a helium- hydrogen composition gradient), the wave excitation mechanism and the 
resultant amplitude of the gravity waves may be quite different from their higher-mass CO counterparts. We compute the 
energy and angular momentum fluxes carried by the waves, and the consequent effect on the orbital and spin evolution of He 
WDs in compact systems. 

We also attempt to calculate the observational signatures of tidal heating in both CO and He WDs. Although our 
calculations of the wave amplitudes are performed in the linear theory, we use these amplitudes to estimate where the waves 
become non-linear as they propagate toward the surface of the WD. If the waves become sufficiently non-linear, they will 
overturn the stratification and break, depositing their energy and angular momentum into the surface layers of the star. Also, 
the WD envelope can develop differential rotation due to the deposition of angular momentum carried by the gravity waves. 
This in turn produces a "critical layer" at which wave absorption takes place due to corotation resonance. We examine the 
criteria for non-linear wave breaking and for wave absorption at the critical layer, and estimate where the outgoing waves 
deposit their energy in He and CO WDs. We then evolve WD models using the MESA stellar evolution code (Paxton et al. 
2011) including a tidal heating term (which is a function of both radius and time) calculated according to various criteria. 
The evolution code allows us to monitor changes in the WD properties as a function of orbital period and thus allows us to 
make predictions of observational signatures for WD binaries as a function of their orbital period. 

The paper is organized as follows. In Section 2, we calculate the amplitude of gravity waves in a He WD model as a 
function of the tidal forcing frequency, and compare these results to those previously obtained for the CO WD models. In 
Section 3, we compute the orbital and spin evolution of a He WD in a compact binary system and estimate the magnitude of 
the tidal heating. In Section 4, we estimate where in the WD the waves will undergo non-linear wave breaking or experience 
wave absorption at a critical layer, and how the location depends on the orbital period and the internal structure of the 
WD. In Section 5, we evolve WD models under the influence of tidal heating and predict observational signatures of the tidal 
heating. Finally, in Section 6, we compare our predictions to observed systems, and we discuss the uncertainties in our results 
and how they may be remedied by future studies. 



2 TIDAL DISSIPATION IN HELIUM WDS 

2.1 Wave Dynamics 

To calculate the amplitude of the tidally excited gravity waves in a He WD model, we use the same method described in 
Paper II. Here we review only the basic concepts and introduce our notations. 
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The Lagrangian displacement associated with the dominant (quadrupole) component of the tidal potential has the form 

€(r,t)= [Ur)r + $±(r)rVx]Y 22 (e i( l>)e- iut , (1) 

where f r is the radial component of the displacement and fx is the perpendicular displacement, r is the radius, Y22 is the 
I — m = 2 spherical harmonic, and uj is the tidal forcing frequency. In this paper, we consider the orbit of the companion 
to be circular and aligned such that ui — 2(f2 — fi s ), where Q is the angular orbital frequency and fl s is the angular spin 
frequency. The displacement £ can be further decomposed into an equilibrium component £ oq that describes the quasi-static 
ellipsoidal distortion of the star, and a dynamical component £ dyn which describes the non-equilibrium wavelike response of 
the star to the tidal forcing. Appendix A describes the details (improving upon the treatment of Paper II) of decomposing 
the equilibrium and dynamical components. 

We calculate the waveform of the tidally excited gravity waves by solving the linear inhomogeneous wave equations for 
stellar oscillations. At the center of the star, we impose the regularity boundary condition. At the outer boundary near the 
stellar surface, we use the (outgoing) radiative boundary condition: 

-(pr 2 /k r )' 
l(^pr 2 /k 



dr 



dyn 



(2) 



1(1 + 1)(N 2 - uj 2 ) 



where p is the density, k r is the radial wave number, 

r „'-■'/■- 

and N is the Brunt- Vaisala frequency. Equation [3] is valid as long as uj 2 <C Lf, where 

t2 _ l(l + l)c 2 



(3) 



(4) 



is the square of the Lamb frequency, and c„ is the sound speed. 

Upon solving the oscillation equations to calculate f r (f) and £,±(r) for a given value of u), the angular momentum flux 
carried by the wave is 

j, (r) = 2muj 2 pr i Rc [i£ yn * £ yn ] . (5) 

When evaluated at the outer boundary r — r out , equation ([5| represents the rate at which the dynamical tide deposits angular 
momentum into the WD envelope. The outgoing angular momentum and energy fluxes can be written as 

Ttidc = J* (r ou t) = ToF(u), (6) 

and 

E tidc = fiToFH, (7) 

where 

GM' 2 ( R s ' 



T 



(8) 



M' is the mass of the companion star, and a is the orbital semi-major axis. The dimensionless function F(u) describes 
the magnitude of wave excitation in the WD, and is strongly dependent on the internal structure of the WD and the tidal 
frequency u. In terms of the commonly used parameterization of tidal dissipation (Goldreich & Soter 1966; Alexander 1973; 
Hut 1981), F(uj) is the related to the tidal phase lag and tidal Q by F(lj) = 3fc25i ag = 3fc2/(2Q) (assuming SI > Q s ), where 
fe is the tidal Love number. 



2.2 Wave Excitation in He WDs 

We perform our calculations on an M = 0.3Mq He WD with a ~ 1.2 x 1O _3 M0 hydrogen envelope, generated using the MESA 
stellar evolution code (Paxton et al. 2011). We evolve the same WD to surface temperatures of T c s = 18000K, T e g = 12000K, 
and T fi = 6000K, with respective radii of R = 2.0 x 10 9 cm, R = 1.6 X 10 9 cm, and R— 1.4 x 10 9 cm. Figure [l] shows propagation 
diagrams for our stellar models. The spike in iV 2 at r ~ 0.8R is due to the He-H composition gradient. 

We have calculated the fluid displacement £(r) as a function of r for tidally excited gravity waves at many values of u. 
Figure [2] shows the wave function £,±(r) and dimensionless angular momentum flux F(lu, r) = j z (r)/To as a function of radius 
for uj — 3 x 10~ 2 , in units where G = M = R = 1. The location of wave excitation can be determined by examining F(uj, r). 
Below the wave excitation region there exists both an ingoing and outgoing wave such that F(u>, r) m 0, while above the 
excitation region there exists only an outgoing wave such that F(ui, r) w constant. It is evident from Figure [5] that the wave 
is excited at the sudden increase in TV 2 associated with the He-H composition gradient. 

Figure [3] displays F(ui) as a function of uj. It is evident that F(uj) is an erratic, non-monotonic function of ui. The reason 
is that the core of the WD behaves as a quasi-cavity containing in-going and out-going waves. At some frequencies, these 
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Figure 1. Propagation diagrams showing p (thick solid line), N 2 (thin solid line), and L| (dashed line) as a function of radius in our 
M = 0.3M Q He core WD model with a ~ 10 _3 M Q hydrogen shell. The three panels arc for WDs with T aB = 18000K (top), T eS = 12000K 
(middle), and T c g = 6000K (bottom). All quantities are plotted in units with G = M = R = 1. 
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Figure 2. Dynamical tide in our He WD model with T a g = 12000K driven by a companion of mass M' = M, with the tidal frequency 
u) = 3.0 X 10 — 2 . Top: The value of F(ui, r) = J z (r)/To (dark solid line) as a function of radius, with J z calculated from equation 151. All 
values are plotted in units of G = M = R = 1. Bottom: The real part of £^ yn (dark solid line) and imaginary part of £ < J yn (dark dashed 
line) as a function of stellar radius. The value of N 2 has been plotted (dashed purple line) in each panel. 



waves constructively interfere at the cavity boundary (the He-H composition gradient) and create a large outgoing wave in 
the WD envelope. At other frequencies, the waves exhibit deconstructive interference, creating a small outgoing wave in the 
envelope (see Paper II for more details). In Paper II we found that on average, F(oj) has the rough scaling F{uj) oc lo for CO 
WDs, whereas we find here that in a rough sense -FX^) oc oj 6 in our He WD models. Moreover, the magnitude of F(ui) of He 
WDs is smaller by nearly five orders of magnitude over the frequency range of interest. 

The smaller values of F(cj) in He WDs stem from the dynamics of gravity wave excitation at a composition gradient. In 
Paper II, we showed that gravity wave excitation due to a composition gradient in CO WDs has a rough scaling 

F(u) » /a 8 , (9) 
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where the notation x indicates that the quantity x should be evaluated in dimensionless units with G = M = R = 1, and / 
is approximately given by 



f 



N 4 g 2 \H a 

-i^)(4Hf)(f)WW* <°™ 
-<ife)(sHS)(S)~WW 4 '» co ™- ™ 

where H is the pressure scale height. The a and b subscripts indicate these quantities are to be evaluated one wavelength 
above the base of the composition gradient, and at the base of the composition gradient, respectively. The middle line of 
equation ( 10 1 is more appropriate for our He WD models, while the last line is more appropriate for our CO WD models. 
Waves are more easily excited at larger values of r a due to the longer lever arm available to be torqued by the tidal potential. 
They are more easily excited at smaller values of Nb because smaller values in Nb produce longer wavelengths which couple 
better with the tidal potential. 

In our He WD models, the value of f a tends to be smaller than in CO WDs because the hydrogen layer is much thicker. 
Furthermore, the value of tends to be larger in He WDs because they are less degenerate than CO WDs (i.e., they have a 
larger entropy and a larger entropy gradient). The difference in stellar structure coupled with the strong dependence on r a and 
Nb causes the value of / to be orders of magnitude smaller in our He WD models than it is in our CO models. Furthermore, 
the slight dependence of r a and N a on u) (due to the fact that these quantities are evaluated one wavelength above the base 
of the composition gradient) creates a slightly steeper scaling F(uj) fa /w 6 in He WDs. 



2.3 Orbital and Rotational Evolution of He WDs 



Having computed the dimensionless tidal torque F(uj), we can calculate the rate at which energy and angular momentum 
are deposited in a He WD in a compact binary. In our analysis, we assume the WD maintains solid body rotation. The spin 
frequency of the WD evolves as 

T F(uj) 



n s = 



(ii) 



where / is the moment of inertia of the WD. The orbital frequency of the WD evolves due to both tidal dissipation and 
gravitational radiation: 



n = 



3T F{lu) 



+ 



Here, \i is the reduced mass of the binary, and tew 



fia" 2tGW 
\a/a\ is the gravitational wave inspiral time given by 



(12) 



tew = 



5c b a 4 
64G 3 MM'Mt 

3.2 x 10 10 sf M ® 



\MM> 



M t 
2M P 



1/3 



n 



0.1s- 



-8/3 



(13) 



At large orbital periods, the second term in equation ( 12 I dominates the orbital and spin evolution such that Q, 3> fi s 



and tidal effects are negligible. However, because of the strong dependence of F(ui) on oj, tidal spin up becomes important at 
short orbital periods. The critical orbital frequency, f2 c , at which tidal spin up becomes important is determined by equating 
Q ~ 3Q/(2tGw) an d ^s- F° r our He WD models, we find [cf. equations (79) and (83) of Paper II] 

5/21 3 / 19 



n c = 



3k m 



5/3 



fGM\ 



10/ M'M 2 / 3 \Rc 2 J 



(GM\ 1/2 



(7.0 x 10"^s _i ) 



fM[Rf 



3/19 



(14) 



where M t = M + M' and k 
Our He WD models have k 



= 1/ (MR 2 ) 
0.094, k = 



R 4 = i?/(10 4 km), and k .17 = k/0.17. 
12000K, Teff = 6000K, respectively. 



Here, M tl = (M + M')/M & , M x = M/M , 
0.13, and k = 0.16 for T cS = 18000K, T cS 
The scaling shown in equation (141 is an important result of our theory]^] The dependence on M and R shows that 
lower mass, larger radius WDs (such as He WDs) should have smaller critical frequencies. So, naively, one would expect tidal 
dissipation to become important at longer orbital periods for low mass WDs. However, as described in Section [2. 2[ the value 



Equations |l4| and |l7| can be derived for the general case Ty^e = Tofui n , as detailed in Appendix [b[ 
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Figure 4. Evolution of the spin frequency Q s in units of the orbital frequency Q as a function of the orbital period. The solid lines 
correspond to our 0.3Af Q WD model with T off = 18000K (black), T off = 12000K (red), and T e g = 6000K (blue). The dashed lines 
correspond to evolutions using F = 3x 10 _3 aj 6 (black), F = 1.5x 10 -2 & 6 (red), and F = 6x 10 _1 (I) 6 (blue) for the WD models of like 
color. The vertical dotted lines denotes the critical orbital period, 27r/f2 c [see equation l |14[ l], for WDs of like color. The black and blue 
curves are indistinguishable from one another. In these evolutions, M' = M and the WDs initially have Q s = 0. 



of / is much smaller in He WDs. These two effects largely offset one another, causing the value of Q c to be similar in He and 
CO WDs. Thus, an He WD in a compact binary will begin to evolve toward synchronization with the orbital frequency at 
roughly the same orbital period as a CO WD, namely, at orbital periods P « 1 hour. 

Figure [4] displays the spin frequency, Q s , as a function of orbital period for an He WD in a compact binary. At periods 
longer than about an hour, the spin frequency increases slower than the orbital frequency, and thus the value of £l 3 /Q remains 
close to zero. However, when the orbital frequency increases to SI « Q c , the value of f2 s /S! begins to increase, i.e., the WD 
becomes more synchronized with the orbit. The WD does not become completely synchronized, however, but retains a nearly 
constant degree of asynchronization such that Q, — fi s ~ Q c (see Section 8.1 of Paper II for more details). 

We can also calculate the amount of heat dissipated in the WD. The heat dissipated in the WD is not equal to the 
tidal energy transfer rate [-Etido, see equation (JT])] because some of the tidal energy flux is stored as rotational kinetic energy. 
Instead, 



At long orbital periods where Q s /Q <C 1, E^ e 
where Q > Q c (see Paper II), 



-Bheat = -Btide(l — fis/fi). 

Etide, with .Etide given by equation |7| 



-Btid 



for 



(15) 

However, at shorter orbital periods 



(16) 



and we find for our He WD model that 



Ehc 



Eud 



19/18 



(1.2 x 10 37 erg 



-In 7/6 r-1/6 , j-17/6 / , Wn,5/6/ \-l/18 
s ) K o'l7/ M l ( M l) {Mtl) 



Q 



0.1 s- 



65/18 



(17) 



Comparison with equation (91) of Paper II reveals that the amount of heat dissipated in an He WD is the same order of 
magnitude as the heat dissipated in a CO WD. At short orbital periods, this heating rate can be much larger than the WD's 
intrinsic luminosity, and it can thus have a substantial impact on the structure, luminosity, and temperature of the WD (see 
Section 4). 
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Figure 5. The tidal energy dissipation rate -B t ide (solid lines) and the tidal heating rate -Ehcat (dashed lines) as a function of orbital 
period for our 0.3M Q WD model with T cff = 18000K (black), T e g = 12000K (red), and T cft = 6000K (blue). Note that at long orbital 
periods, Etidc — -^heat an d these curves overlap. The dotted back line is the energy flux carried away by gravitational waves, i?GW- I n 
these evolutions, M' = M and the WDs initially have fl s = 0. 



Figure [5] shows the tidal energy flux and heating rate as a function of orbital period for a He WD in a binary with 
M' = M. It is clear that Egw 3> -Etide at all orbital periods, where £qw = G M M' / (2atGw) is the energy flux carried away 
by gravitational waves. At long orbital periods where Q, < Q c , the value of Shcat — -Btide, but the magnitude of the heating is 
small and has a negligible effect on the WD. At shorter orbital periods where Q > Q c , the heating rate is well described by 
equation ( 17 1. 



3 LOCATION OF TIDAL HEAT DEPOSITION 

Until now, we have not examined the processes that convert the angular momentum and energy carried by gravity waves 
into the rotational angular momentum and internal energy of the stellar envelope. We have assumed that the gravity waves 
propagate into the envelope of the WD where they somehow dissipate, depositing their energy and angular momentum. Here, 
we investigate these processes so that we can estimate the magnitude of tidal heating as a function of depth within the star. 



3.1 Non-linear Wave Breaking 

As the tidally excited gravity waves propagate outward in the WD envelope, their amplitudes increase. In the WKB limit, with 
u) <5cC iV and uj Li, the radial wave number is given by k r ~ —k±N/uj, where the horizontal wave number is kj_ = y/UJ, + l)/r. 
The amplitudes of the radial and horizontal displacements scale as 

c ikrr t N tia\ 

Obviously, |fc r |/fc_L = N/u) > 1 and |£ r /£±| = y/l{l + Vju/N < 1. 

As the gravity waves reach sufficiently large amplitudes, they are expected to break and quickly damp, locally depositing 
their energy and angular momentum. The critical amplitude for wave breaking may be estimated by comparing the Eulerian 
acceleration dv/dt = — cj 2 £ (where v = — iu)£ is the fluid velocity) with the non-linear "advective" term v • Vv = — lj 2 £ ■ V£. 
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For gravity waves, most of the terms in the advective derivative satisfy the non-linearity condition £ ■ V£ ~ £ when 



1(1 + 1) 



6 



l. 



(19) 



Thus, nonlinear effects become important when the radial (horizontal) displacement is comparable to the radial (horizontal) 
wavelength. Alternatively, one may expect the waves to begin breaking when the shear is large enough to overturn the 
stratification of the star, i.e., when the Richardson stability criterion, N 2 /\dv±/dr\ 2 > 1/4, is violated. This occurs when 
i5 1, a condition similar to equation ( f9 1 for I — 2. Equation ( 19 1 is similar to the wave breaking criterion discussed in 
Ogilvie & Lin (2007) and found in three-dimensional simulations of gravity waves approaching the center of a solar-type star 
(Barker & Ogilvie 2011). 

However, equation (19 1 corresponds to a physical displacement, [£| ~ |£x| ~ R/[l(l + 1)] in the envelope of the star, ft 
seems unlikely that fluid displacements of order the radius of the star can be realized before non-linear wave breaking occurs. 
Moreover, the f component of the advective derivative contains the term — £ ± ■ $j_/r. This implies that non-linear effects may 
become important when £±/r ~ £ r , which corresponds to fc r £± ~ 1(1 + 1) in the WKB limit. Thus, we will also consider the 
non-linear breaking criterion 



(20) 



where /3 <; 1 is a free parameter. Choosing different values of j3 will allow us to test how the location of tidal energy deposition 
depends on different wave breaking criteria. 

To implement the non-linear criterion (1191) or |20l for different orbital frequencies and companion masses, we can use 



our numerically computed wave function £(r) and extend it to the near surface region via the WKB amplitude relation ( 18 I 
Alternatively, we can use equations jHp, and (18 1 to find 

F(u) Q. 4 N(r 1 ' 1 



Mx.(r)\^R 



(21) 



2m^/l(l + 1) Gp(r)u; 

with I = m = 2. The radial displacement can be obtained from equation ( |18[ |. The above expression gives the wave amplitude 
in the outgoing wave region, i.e., the region above the composition gradient at which the waves are excited. 

We denote the location of non-linear wave breaking by AMb, the mass above the point at which the wave amplitude 
satisfies one of the non-linear criteria discussed above. Figure [6] shows AMb as a function of orbital period for both CO and 
He WD models, calculated according to equations ( 19 I and ( |20[ ). At large orbital periods (P > 40 minutes), the waves do not 
become non-linear under criterion (191 below the surface convection zone (whose depth is AM conv — 1O _12 M0 for the CO 
WD and AM conv ~ 1O _13 M0 for the He WD). In this case, the g-mode analysis of Paper I and Burkart et al. (2012) may 
become applicable. At shorter orbital periods, the waves become non-linear according to equation (19 1, but they break near 
the surface of the WD where AM B < 10" 8 M Q . 

However, under the criterion of equation (20) with f) = 1, the waves become non-linear deeper in the star, and nearly 
always reach non-linear amplitudes before encountering the surface convection zone. In the CO WD model, the value of AMb 
jumps upward to AMb ~ 1O _4 M0 at an orbital period of about 30 minutes, while AMb stays below 1O _5 M0 in the He 
WD. The reason is that the CO WD has two composition gradients whereas the He WD has only one. The He-H composition 
gradient in the CO WD causes the outgoing waves excited at the C-He gradient to have larger values of both |£x| and \k r \ 
in this layer, promoting non-linearity via equation (20 1. No such composition gradient exists above the excitation region in 
an He WD, so the waves do not reach non-linear amplitudes until they are close to the surface. For reference, Figure [6] also 
shows the value of AMb calculated with the intermediate non-linearity conditions |fc r £j_| = 3 and |fc r £x| = 10. Using these 
more conservative criteria yields results similar to the criterion |fc r £x| = 1, but with generally smaller values of AMb- 



3.2 Wave Absorption at a Critical Layer 

Until now, we have considered the WD to be rotating as a rigid body. Our results have indicated that in a rigidly rotating 
WD waves will break and deposit their angular momentum in the outer layers of the star where AMb < 1O _4 M0. Since a 
small fraction of the stellar mass absorbs the entirety of the angular momentum flux carried by waves, the outer layers of the 
star may spin up rapidly. If the outer layer spins up faster than angular momentum can be transported to the core, it will 
attain synchronous rotation with the orbit of the companion. The outgoing gravity waves will then encounter a critical layer 
(corotation resonance), where the wave frequency in the rest frame of the fluid, to = 2[fi — f2 s (r)], equals zero. At the critical 
layer, the gravity waves are directly absorbed by the fluid since their radial wavelength |fc r j _1 ~ to/(k±N) and group velocity 
v g ~ uj 2 /(k±N) both approach zero. As the continuous train of gravity waves deposit their angular momenta in the stellar 
envelope, the critical layer will move to larger depths. Such an "outside-in" scenario for tidal spin-up was first discussed by 
Goldreich & Nicholson (1989) for massive (early- type) stellar binaries, and was applied to WD binaries in Paper III. 

In this section, we model the spin evolution of the WD using a simple two-zone model. In this model, the envelope 
of the star forms a critical layer that rotates synchronously with the orbit (f2 Sj0n v = fi), while the core of the star rotates 
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Orbital Period (minutes) 



Figure 6. The envelope mass AMj above the point at which outgoing waves become non-linear as a function of orbital period for an 
M = O.6M CO WD model with T cff = 10000K (top) orbiting an M = O.3M He WD model with T cff = 12000K (bottom). The solid 
black lines are calculated from equation \19\ , while the three dashed lines are from equ ation ||20[ l with /3 = 1, 3, 10 (from top to bottom). 
This plot assumes the WD spin and orbital frequencies evolve according to equations jll|12|. 



sub-synchronously (f2 s ,corc < fi). The envelope and core are coupled, with angular momentum being transferred to the core 
according to a parameterized coupling time, t coup . The angular momentum of the system evolves according to 

Jcnv — ^(^corc) ~ (^s,env ^s,corc) (^^) 

tcoup 

'/core — 7 (^s,cnv ^s.corc)- (^^) 

Ccoup 

Here, j z is the angular momentum flux which can be calculated from equation (Jfjl. We assume that the waves are excited in 
the core and dissipated in the envelope. Consequently, the angular momentum source term J z is only present in the envelope 
evolution equation, although it is dependent on the tidal frequency in the core, c^Wo = 2(f2 — f2 s , corc ). Using f2 SlCn v = Q and 
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the fact that the star's total moment of inertia is / = / cnv + /core, equations (22 1 and (23) may be rewritten 



SI 



SI -SI. 

CO 



I-Ie, 



^coup 



(24) 



(25) 



Since the orbital decay is dominated by the emission of gravitational waves, Sl/Sl ~ ?>Sl/(2tc 



With appropriate initial conditions, we can integrate equations (24 1 and (251 to calculate the values of 7 onv and fi SiCore as 
a function of orbital period. We then obtain the mass AM cnv of the envelope corresponding to the value of / on v In this simple 
two-zone model, the tidal heat is deposited entirely at the base of the envelope where AM = AM cnv . The thickness of the 
envelope is dependent on the parameter t coup . Unfortunately, angular momentum transport in stars is not well understood. 
In stably stratified stars like WDs, angular momentum is likely transported by magnetic torques, e.g., via the Tayler-Spruit 
dynamo (Spruit 2002). In Appendix |c| we estimate value of t, 
parameters. 



coup, and find that lyr < t coup < 10 4 yr for realistic WD 
We can calculate approximate equilibrium solutions to equations (24 1 and (25 I. Since t coup <C tGw, the first term on the 



right hand side of equation (241 is negligible. Then at equilibrium, 



core / L c 

si-si s 



(26) 



At large orbital periods where SI <Sl c [see equation ( 14 1] the system is not synchronized such that Sl 3: 

<Jz(2Sl)t C ou P 



-C SI and 



SI 



for 



At short orbital periods (SI > Sl c ) where synchronization has begun, SI 
3ISl/(2tavf), and we have 



SI < Sl c 

Sis. 



2f2 c tGw 



for SI > Sl c 



(27) 

Sl c . In this regime, J z = E t ide/Sl ~ 

(28) 



Comparison with our numerical integration of equations ( 24 1 and ( 25 \ shows that the approximations of equations ( 27 1 and 



(28 1 are very accurate. 

Figure [7] plots the value of AM H1V as a function of orbital period for our two WD models, using values of t coup ranging 

= 0, as is appropriate at long 



and Sl B 



the value of AM env remains small at all orbital 



from lyr to 10 yr. We begin our calculation at P or b > lhr and use 7 el 
orbital periods where tidal effects are small. For the chosen values of t cc 
periods (AA/ cnv < 10 _2 Mq). Thus, if a critical layer develops in a real WD, we expect it to be restricted to the outer region 
of the star. However, the critical layer extends to very large optical depths, suggesting that binary WDs may be observed to 
be synchronized at large orbital periods even if their cores are not synchronized. Our results indicate that, for our CO WD 
model, the critical layer most likely does not penetrate as deeply as the C-He composition gradient where gravity waves are 
excited (or the He-H gradient in our 0.3Mq He WD model), so our assumption that J z is a function of o; corc = 2(Sl — Sl s ,coic) 
is justified. Finally, we note that the values of AM cnv are similar in magnitude to the values of AMb calculated in Section 

ED 



4 THERMAL EVOLUTION OF TIDALLY HEATED BINARY WHITE DWARFS 

In Section [2] (and in Paper II), we demonstrated that the tidal heating rate -Ehcat of a WD in a compact binary may be 
substantially larger than the intrinsic luminosity of the WD. However, the consequence of this tidal heating is not clear a 
priori and depends on the location of heat deposition (Section |3|. The assumption that the WD reaches a state of thermal 
equilibrium is not justified because the WD cooling time can be much larger than the gravitational inspiral time. In this 
section, we calculate the effect of the tidal heating on the WD structure and on its surface temperature and luminosity. 



4.1 Details of Calculation 

To calculate the effect of tidal heating, we first calculate the value of AMb according to the wave breaking criteria of Sections 



3.1 and 3.2 (for the purposes of this calculation, we refer to the value of AM mv as AMb, because it determines the depth 



at which tidal heat is deposited). We deposit the tidal heat uniformly per unit mass in the outer layers of the WD that have 
AM < AMb- Although the radial dependence of this heating function is unlikely to be realistic, we find that the results are 
not strongly dependent on the form of the radial heat deposition (although they are sensitive to the value of AMb)- The 
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Figure 7. The envelope mass AA/ env above the critical layer as a function of orbital period for our M = 0.6Mq CO WD model with 
T cff = 10000K (top) orbiting our M = O.3M He WD model with T cff = 12000K (bottom). The solid black lines have i coup = lyr, 
the dot-dot-dashed green lines have t coup = lOyr, the dot-dashed orange lines have t coup = 10 2 yr, and the dashed red lines have 
tcoup = 10 3 yr. 



heating rate per gram of material, £hcat, is then 

e hcat = for AM > AM B (29) 

ehcat = for AM < AM B , (30) 

IXlvl B 



with the value of iShcat calculated from equation ( 171 (and its counterpart in Paper II for CO WDs). 

To understand the effect of tidal heating on the WD properties, we evolve WD models using the extra heating term 
calculated via equation < |30[ ). We use the one-dimensional stellar evolution code MESA (Paxton et al. 2011) to evolve our 
WD models. At the beginning of the evolutions, the WDs have the same profiles used to calculate the magnitude of wave 
excitation (e.g., the models shown in Figure [T] for our He WD models). The initial orbital period is one hour, with companion 
masses discussed below. 

During the course of our evolutions, we do not calculate new values of F(oj) and fi c at each time step. In principle, these 
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values change as the stellar structure adjusts to tidal heating (and WD cooling). However, since waves are excited at a depth 
well below where they deposit their energy as heat, the properties of the WD at the location of wave excitation experience 
little change during the evolution. Hence, we expect the values of F(uj) and Q c to remain roughly constant. In contrast, the 
value of AMb is strongly dependent on orbital period, and is updated at each time step. 

In general, the amount of tidal heat depends on the masses of the two WDs. We do not attempt to cover the whole 
spectrum of WD masses. Instead, we consider only two cases: a 0.6Mq CO WD paired with a O.SMq He WD (in this case we 
examine the heating of both WDs), and a O.6M CO WD paired with a 0.9M Q companion star (in this case we examine only 
the O.6M WD). We present results for the 0.6M @ CO WDs at initial temperatures of 5000K, 10000K, and 15000K, and for 
the 0.3Mq He WDs at initial temperatures of 6000K, 12000K, and 18000K. These temperatures roughly span the observed 
temperatures in compact WD binaries (Kilic et al. 2012). 



4.2 Effects of Tidal Heating 

Figure[8]displays the surface temperature of the 0.3Mq He WD as a function of orbital period. For comparison, we also plot the 
temperature of a non-tidally heated WD (i.e., we set i5 nC at = in our evolutions). At long orbital periods (P > 30 minutes), 
the tidal heating has little effect on the luminosity or temperature of the WD, except to slightly raise the temperature of the 
To = 5000K model. In this phase of its evolution, the WD remains well described by the cooling track for a non-tidally heated 
WD. At shorter periods (P < 20 minutes), the luminosity and temperature become substantially larger due to tidal heating. 
At the smallest orbital periods [P < 8 minutes), the WD has a temperature in excess of 18000K regardless of its initial 
temperature, and the luminosity is dominated by escaping tidal heat. The case with t coup — 10 yr and initial temperature 
Teff = 6000K (top right panel) ends abruptly, indicating that the surface hydrogen layer ignited to create a tidally induced 
nova, which we discuss in greater detail below (see also Paper III). 

The top left panel of Figures [8] [9] and [To] are calculated for cool WDs with deep convection zones. In these WDs, the 



waves do not become non-linear below the convection zone according to equation ( 19 1 at any orbital period, and they may 
become non-linear according to equation ( |20[ ) at orbital periods below one hour. When the waves do not become non-linear, 
we distribute tidal heat evenly throughout the WD, but a more accurate calculation should use the discrete g-mode formalism 
discussed in Paper I (c.f. Burkart et al. 2012, Valsecchi et al. 2012). 
Figures |8|10| also show the tidal heating temperature, defined as 

_ / -Bhe ' 



TcS ^ = [^k) ' (31) 

where a is the Stefan-Boltzmann constant. In the 0.3M© WD, the tidal heat diffuses to the surface very quickly (see Section 



4.3 for more discussion), regardless of the non- linear breaking criterion. In this case, the observed luminosity is roughly 
L ~ Li n t + Sheat, where Li nt is the intrinsic luminosity of the cooling, non-tidally heated WD. 

Figures [9] 10| display the temperature of a 0.6Mq WD as a function of orbital period, for companion masses of 0.3Mq and 
0.9Mq. For the cases using the non-linear breaking criterion of equation ( 19 1 and critical layer absorption with t C ou P = lyr, the 



results are very similar to those of the O.3M0 He WD models. When the luminosity is dominated by tidal heating, the 0.6Mq 
WD models have slightly larger surface temperatures than the O.3M0 models. Although the value of -EWat is comparable 
between the O.3M0 and 0.6Mq models, the smaller radius of the 0.6Mq WD models requires a larger surface temperature to 
radiate the same amount of energy. Finally, the luminosities and temperatures are usually larger for the 0.9M Q companion, 



as expected from the scaling of equation (17 1 



However, for the non-linear breaking criterion of equation (201 and the critical layer absorption with t C ou P = 10 3 yr, the 



results are markedly different. For these criteria, most of the tidal heat is deposited deeper in the WD where it cannot quickly 
diffuse outward to be radiated away. These criteria thus lead to generally lower WD surface temperatures. However, because 
the tidal heat is not quickly radiated away, the layers in which the heat is deposited may heat up substantially. 

The layers just above the He-H composition gradient are primarily composed of degenerate hydrogen. If these layers are 
able to trap enough heat, their temperature will increase until hydrogen fusion begins. Due to the degeneracy of the hydrogen, 
the ignition of fusion starts a thermonuclear runaway similar to a classical nova. Our evolutions show that these tidally induced 
novae occur in our two cooler O.6M0 CO WD models and in our coolest O.3M0 He WD model. The novae occur only if the 
gravity waves deposit their heat near the base of the hydrogen layer, i.e., only for the heat deposition criteria which yield 
10 _5 M© < AM B < 10~ 3 M Q for our CO WD model and 10 _4 M© < AM B < 10- 2 A'/ Q for our He WD model. Novae do not 
occur for the warmer models because the hydrogen is not degenerate, so these models are able to burn the hydrogen stably. 

In a thermonuclear runaway event, most of the hydrogen will be burned to helium or will be ejected from the system 
(Truran 2002). We do not attempt to predict a detailed observational signal of such an event, other than to speculate that 
it will be very similar to a classical nova. The thermonuclear runaway may dramatically change the dynamics of subsequent 
tidal heat deposition. Assuming the hydrogen shell is much thinner after the nova event, tidal heat may be deposited closer 
to the surface where it can quickly diffuse outwards, similar to the results for our warm He WDs. 
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Figure 8. The surface temperature of our 0.3Mq He WD model with a 0.6A/q companion as a function of orbital period, for initial 
temperatures of 6000K (top), 12000K (middle), and 18000K (bottom). On the left-hand side, the solid black lines are calculated with 
the wave breaking criterion of equation (19) , while the dashed red lines are calculated with equation ([20J . On the right-hand side, the 
solid black lines are calculated via the two-zone model with t C oup = lyr, while the dashed red lines are calculated with t coup = 10 3 yr. 
The dotted lines are calculated for a WD with no tidal heating and the same initial temperature. The blue dot-dashed lines are the 
values of and T e g ^ c from equation (31| . The black and red lines overlap at large orbital periods, while the black and blue lines often 
overlap at small orbital periods. Discontinuities in temperature are due to sudden changes in the location of heat deposition (see text). 
Stars indicate the occurrence of a tidal nova. The plot extends to an orbital period of 5 minutes, the approximate orbital period at which 
Roche lobe overflow occurs. 



4.3 Heating and Cooling Time Scales 

The effect of tidal heating can be better understood by examining the relative time scales of the WD inspiral, heating, and 
cooling processes. The WD inspiral time due to gravitational radiation, £gw, is given by equation (131. The inspiral time scale 
ranges from more than a Hubble time at large orbital periods (several hours), to less than a million years at short orbital 
periods (less than about fifteen minutes). Any process that acts on a time scale longer than the inspiral time is irrelevant, 
because the WDs will have merged (or begun stable mass transfer) by the time the process makes a substantial change. 

The temperature evolution of the WD is determined not only by the amount of tidal heat deposited, but also by the rate 
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Figure 9. Same as Figure [8] but for our 0.6M© CO WD model with a 0.3A/q companion. 



at which that heat is able to diffuse to the surface. The time scale on which a shell of material heats up (in the absence of 
cooling) is 



c T 

ihcat(V, fi) = — , 



(32) 



where c p is the specific heat at constant pressure. The radial profile of the heating time depends primarily on the magnitude 
and radial dependence of the tidal heating. As tidal heat is deposited in a shell of the WD, it will diffuse on a thermal time 
scale, 

Pc n T 



£therm(r, fi) = 



(33) 



where P is the pressure and F is the heat flux through the shell. The thermal time scale has a very sensitive dependence 
on the depth of the shell in question (see Figure 111. In the core of the WD, itherm ~ 10 9 years, whereas near the surface, 
ithorm < 1 year. 

Figure [TT] shows a plot of Jgw, tteati and ttherm as a function of AAf for our Q.QMq WD with a 0.3Mq companion at 
orbital periods of 45 and 15 minutes. At long orbital periods, tthcrm <fC ihcat at all radii. The tidal heat is able to quickly 
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Figure 10. Same as Figure[8] but for our 0.6Mq CO WD model with a 0.9Mq companion. 



diffuse to the surface and be radiated away. The temperature of the WD reaches a thermal equilibrium such that tidal heat is 
radiated at the same rate it is deposited. We find this is also the case for our 0.3Mq He WD model at most orbital periods. 

At short orbital periods, ttherm ~ theat near the base of the heat deposition zone. The temperature profile will adjust 
so as to re-establish thermal equilibrium such that ttherm < theat, thereby changing the internal structure of the star. When 



AMb is calculated with equation (191, the star is able to adjust to the heating by steepening its temperature gradient (thus 
increasing its luminosity) such that ttherm < thcat at all radii [^However, when AMb is calculated with equation (20 1, the star 
is unable to reach thermal equilibrium and ttherm > thcat at the base of the heat deposition zone. The shell at AM w 1(T 4 M 
heats up, and eventually it reaches temperatures high enough to trigger run-away hydrogen fusion. 

In principle, tidal heating may change the structure of the WD enough to alter the dynamics of wave excitation and wave 



2 In Figure [TT] the value of ttherm goes to infinity at some values of AM because the heat flux goes to zero at these locations. This can 
occur when the amount of tidal heat diffusing inwards from the surface is equal to the amount of intrinsic WD heat diffusing outwards 
from the core (i.e., there is a local temperature minimum). It can also occur where the tidal heat diffuses equally in both directions (i.e., 
there is a local temperature maximum). 
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Figure 11. The thermal time [solid lines, equation j33l ], heating time [dashed lines, equation ( |32| ] , and inspiral time [dotted lines, 
equation | |13[ )] as a function of exterior mass, AM, in our 0.6A/q CO WD model with initial temperature T c ff = 10000K. The black lines 
are calculated using the breaking criterion of equation l |19[ ), while the red lines are calculated using the breaking criterion of equation 
\20) . The top panel is calculated for an orbital period of 45 minutes, while the bottom is for an orbital period of 15 minutes. Both panels 
are calculated for a 0.3Mq companion. 



breaking. However, we find that this is not the case, with the exception of the instance of a run-away hydrogen fusion event. At 
and below the location of wave excitation (AM > 6 x 10 -3 ), tow <IC theatj so these layers are not affected by tidal heating. At 
shallower depths (AM < 6 x 10 -5 Mq), significant heating may occur, creating large temperature gradients. However, in our 
evolutions no interior convection zone forms, despite the large temperature gradients. We thus conclude that the locations of 
wave excitation and breaking will not be significantly altered by the tidal heating, except in the case of a run-away hydrogen 
fusion event. 



5 DISCUSSION 

We have calculated the amplitude of tidally excited gravity waves in a 0.3Mq He WD in a compact binary system, using an 
outgoing radiative boundary condition. This amplitude translates directly into the rate at which tidal energy and angular 
momentum are transferred into the WD as a function of the tidal forcing frequency. As in the case of CO WDs (see Paper 
II), we find that the tidal torque and energy flux depends on the tidal forcing frequency in an erratic way. On average, the 
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dimensionless tidal dissipation rate F(w) [related to the inverse of the tidal quality factor, see equation jfjj], is several orders 
of magnitude smaller for He WDs than for the O.GMq CO WDs studied in Paper II. This difference arises from the larger 
entropy gradients in low-mass WDs, causing decreased coupling between gravity waves and the tidal potential [see equation 
( 10 1] . Nevertheless, since low-mass He WDs have larger radii, we find that the critical orbital frequency above which spin-orbit 
synchronization starts is similar in He and CO WDs, both occurring at an orbital period of about an hour. Furthermore, 
the amount of tidal heat deposited is similar in He and CO WDs, exceeding the intrinsic luminosity at short orbital periods 
(P < 30 minutes). 

We have also estimated the location of tidal heating due to the non-linear breaking of outgoing gravity waves or resonant 
absorption at a critical layer in the WD envelope. The tidal heat is likely deposited in the outer layers of the WD with 
AM < 1O" 2 M , although the precise location depends on the details of non-linear wave breaking and the efficiency of 
rotational coupling between the WD core and surface layers. We have evolved WD models to evaluate the effect of tidal 
heating as function of orbital period. Tidal heating is unimportant at long orbital periods (P or b > 1 hour), but can dominate 
the luminosity of the WD at small orbital periods (P or b <■ 15 minutes). We have also found that if the tidal heat is deposited 
deep enough in a WD (near the He-H composition gradient), it may build up enough to trigger a nova-like hydrogen burning 
event (a "tidal nova"). 

Our theory can be constrained by comparing the predictions of our tidal heating calculations to observed compact WD 
binaries. The two 39 minute systems J0106-1000 (Kilic et al. 2011a) and J1630+4233 (Kilic et al. 2011b), and the 12.75 minute 
system J0651+2844 (Brown et al. 2011) provide the best opportunities. We find that the warm temperatures of the He WD 
primaries in the two 39 minute systems (T ff = 16490K in J0106 and T c g — 14670K in J1630) can not be explained by tidal 
heating. These WDs are likely young or kept warm by residual hydrogen shell burning (Steinfadt et al. 2010b). The unseen 
companions in these two systems do not have measured effective temperatures, but they are unlikely to be substantially 
increased by tidal heating. 

The 12.75 minute system J0651 likely exhibits strong tidal heating signatures. This system is composed of a primary He 
WD with T cff = 16530K and M = O.26Af , and a secondary CO WD with T eff w 8700K and M = 0.50M© (Hermes et al. 
2012). Comparison with Figure [8] indicates that the luminosity of the primary is mostly due to intrinsic heat, but that tidal 
heat may slightly increase its temperature and luminosity. In contrast, comparison with Figure [12] indicates the luminosity 
of the secondary is likely dominated by tidal heat. In fact, Figure [P2| predicts that T e g is typically larger than the observed 
surface temperature, except for relatively deep tidal heat deposition. 

Our results indicate that the secondary in J0651 is fairly old (it takes a WD several 10 9 yr to cool to T e ft = 5000K), and 
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that tidal heat deposition occurs relatively deep in the star, at AMb > 10 _4 Mq. Such deep tidal heat deposition could arise 
due to a thick critical layer created by inefficient rotational coupling between the surface layer and the core, or it could be 
due to non-linear wave breaking at the He-H composition gradient. Furthermore, Figure [f2] indicates that a tidal nova may 
occur in the future of J0651. 

However, there are other possible explanations for the cool observed surface temperature of the secondary in J0651. One 
possibility is that the value of -Eheat used in our calculations is too large. This would require the value of Q, c to be smaller than 
suggested by our calculations, i.e., tidal effects become important at longer orbital periods than we have predicted. Stronger 
tidal effects may be possible if some other tidal dissipation mechanism (e.g., inertial waves, non-linear wave excitation, or 
spin-up via resonance locking, see Burkart et al. 2012 and discussion below) plays an important role in tidal synchronization. 
This would cause the tidal heating rate at shorter periods to become smaller, which could explain the low observed value of 

Teff- 

We find that the ignition of run-away hydrogen fusion is a general characteristic of the tidal heating process for CO WDs 
with hydrogen envelopes, provided that a significant fraction of the tidal heat is deposited near the base of the hydrogen 
shell. Thermonuclear runaway only occurs in WDs with initial surface temperatures T c g < 1.2 x 10 4 K, otherwise tidal heating 
promotes steady hydrogen burning. The thermonuclear runway usually occurs at orbital periods 5min < P or b < 20min, 
depending on the location of heat deposition, initial temperature of the WD, and companion mass. We speculate that the 
thermonuclear runaway will create an event very similar to a classical nova. Since most of the hydrogen would likely be burned 
or ejected in such an event, recurrent novae are unlikely. Thus, the occurrence rate of these tidally induced novae may be 
comparable to that of WD mergers involving a CO WD. 

Finally, we note that we have ignored the effect of mixing in our WD evolutions. Substantial mixing may be caused by 
the turbulence of breaking gravity waves or by shear instabilities due to differential rotation. If mixing occurs faster than the 
gravitational settling time, the composition gradients inside the WD may be smoothed out. Since the composition gradients 
play important roles in the excitation and breaking of gravity waves, substantial mixing may change the dynamics of the tidal 
synchronization and heating process. We have also ignored the effect of crystallization on the wave dynamics. Although this 
issue is unlikely to affect He WDs because of their long cooling times, it may be important for cool CO WDs. 

Soon after our paper was submitted, Burkart et al. (2012) posted an article on arXiv studying tides in WD binaries. In 
addition to examining gravity waves that break [using the criterion of equation ( 19 1] in systems with short orbital periods, they 
discussed resonance locking with WD g-modes at larger orbital periods (P > lhr). This may cause substantial tidal spin-up 
of the WD in the orbital period range lhr < P < 4hr if the associated gravity waves do not break. At shorter orbital periods 
(P < 45min), Burkart et al. (2012) found that gravity waves indeed break, producing tidal torques and spin evolution similar 
to ours. Note that if gravity waves break via the non-linearity criterion of equation (201, then the resonance locking regime 
will be limited to long orbital periods. This will reduce the effectiveness of tidal spin-up prior to the onset of continual gravity 
wave breaking, bringing their results more in line with ours for P ~ hrs. In another recent paper on WD binaries, Valsecchi 
et al. (2012) discussed possible "anti-resonance locking" of g-modes in the linear theory. The nature of this "anti-resonance" 
locking is not clear, nor is the parameter regime where such locking is effective. Overall, it appears that there is a general 
agreement that at short orbital periods (P < lhr), the WD spins up on the orbital decay timescale, with Q s /£l approaching 
unity as increases, but maintaining an approximately constant (Q — f2 s ) until binary merger or the onset of mass transfer 
[see Figure 4 or equation (B3l; see also Section 8.1 of Paper II]. 

Future observations and simulations can further constrain these theories. We are hopeful that more compact WD systems 
with P or b < 30 minutes will be discovered in the near future. Measurements of the masses, luminosities, temperatures, and 
spins of these systems will provide more data points for comparison with our theory. Finally, simulations of gravity waves 
propagating through a WD envelope could be used to understand the non-linear criterion that governs wave breaking in WDs 
and could more conclusively determine where tidal heat is deposited in a WD. 
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APPENDIX A: CALCULATION OF EQUILIBRIUM TIDE 



In this appendix we present our improved method for isolating the dynamical and equilibrium components of the tidal 
perturbation. Using the Cowling approximation, the oscillation equations are 
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5P' + \SP + p(N 2 - ui 2 )£ r + pU' = 0, 
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where U is the tidal potential produced by the companion, the ' denotes d/dr, and g is the gravitational acceleration. The 
other perturbation variables are related to SP and £ r by 
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Using equation ( A3 1 , equations (|A1| and (|A2j may be rewritten 
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The zeroth order solution to the equilibrium tide can be found by taking the limit lo = in equations ( A5 \ and ( A6 1 
yielding (see also Goldreich & Nicholson 1979) 
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Defining f r = fj: + £ r and likewise for fx, we substitute into equations |A5| and (A6l to find 
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To find the first order term of the equilibrium tide, we again take the limit of to = in equations ( A9|A10 l to find 
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We find that higher order terms are not essential for our purposes. Moreover, these terms typically involve multiple 
derivatives of stellar properties that are difficult to compute from a grid of stellar quantities. The total equilibrium tide is 
£ Cq = £ Cq ,o + £ Cq ,i + (u A ). M our outer boundary, the dynamical tide is computed using f dyn = £ - f cq . 



APPENDIX B: GENERAL SCALING FOR SPIN-ORBIT SYNCHRONIZATION AND TIDAL HEATING 

In Section 2.3, we showed that the averaged dimensionless tidal torque on a He WD scales with the tidal forcing frequency as 
F oc oj 6 . For a CO WD, the scaling is F oc lu 5 (see Paper If). In general, we may parameterize the tidal torque by 

T tidc = j, = T fiu n (Bl) 

with the dimensionless values / and n determined from the dynamics of gravity wave excitation and dissipation within the 
WD. Using equation (Bl I, we can examine the general behaviors of spin-orbit synchronization and tidal heating in a compact 
binary undergoing gravitational radiation-driven orbital decay (see also Section 8.1 of Paper II). 
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Combining equations (|6| and ( 12 I as we did in Section 2.3 we find that tides begin to synchronize the star at a critical 
orbital frequency 
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Here, R s = 2GM/c 2 , q M = (M + M') 5/s /{M 2/3 M'), and fi dyn = ^GMjR?. Thus, less compact stars (with larger i?/i? s ) 
begin synchronization at larger orbital periods. Note that if n <; 3 as for WDs, the critical frequency fi c does not depend 
sensitively on / or k, and only a very rough estimate of these quantities is necessary. 



The tidal heating rate is given by equation ( 15 1. For D. > D. c (see Paper II), 



and combining equations \15) , ( |B2[ ), and ( |B3[ ) yields 
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where g M 2 = (A/ + M') (5+n)/3 M 2(n " 1)/3 /M' n+1 . Note that for n £ 3 as it is for WDs, the tidal heating rate is very insensitive 
to the value of /. At typical orbital periods (SI <C fidyn), the tidal heating rate is a small fraction of the energy loss rate due 
to gravitational waves (especially for non-compact stars), but it becomes more significant at shorter orbital periods. 



APPENDIX C: ESTIMATING THE CORE-ENVELOPE COUPLING TIME SCALE 



In the two zone model (see Section 3.2 1, the thickness of the synchronized envelope is dependent on the parameter t coup . In 
stably stratified stars like WDs, angular momentum can be transported by magnetic fields. In the presence of a poloidal field 
B connecting the core and envelope, t coup can be estimated from the Alfven wave crossing time, 

t A = f\^B. (ci) 

We find t A w 10 2 yr(10 3 G/B) for our CO WD model, and t A ~ 50 yr(10 3 G/B) for our He WD model. 

For WDs without an intrinsic magnetic field, angular momentum may be transported via the Tayler-Spruit dynamo 
(Spruit 2002). To estimate t C ou P , we calculate the effective viscosity vts for angular momentum transport using the the 
method outlined in Spruit (2002). For simplicity, we calculate uts without including the effects of composition gradients in 
the WD to find 

a \ 1 ' 2 f * \ 1 /' 2 



where kt is the heat diffusivity. This viscosity does not depend on the local shear, and is thus independent of the precise 
rotation profile of the star|^] 

We then estimate the coupling time from the Tayler-Spruit dynamo to be 

Its = j R dr?—^. (C3) 

JO V TS 

We find t TS w 2 x 10 3 yr (P/30min) 3/2 for both O.3M He and 0.6M Q CO WDs. The values of t A and t TS convince us that 
icoup < 10 3 yr for the short orbital periods of interest. 



REFERENCES 

Alexander, M.E. 1973, Astrophys. Space Sci., 23, 459 
Barker, A., Ogilvic, G. 2011, |arXiv:1102.0861| 
Bloom, J., et al., 2012, ApJL, 744, 17 

Brown, R.B., Kilic, M., Hermes, J.J., Allende Prieto, C, Kenyon, S.J., Winget, D.E. 2011, ApJ, 737, 23 

Burkart, J., Quataert, E., Arras, P., Weinberg, N., 2012, |arXiv:1211. 1393] 

Dan, M., Rosswog, S., Guillochon, J., Ramirez-Ruiz, E. 2012, MNRAS, 422, 2417 

Di Stcfano, R. 2010, ApJ, 719, 474 

Fuller, J., Lai, D. 2011, MNRAS, 412, 1331 (Paper I) 

Fuller, J., Lai, D., 2012, MNRAS, 421, 426 (Paper II) 



3 The Tayler-Spruit dynamo requires a minimum local shear to operate (see Spruit 2002), causing significant shear to persist if the 
dynamo dominates angular momentum transport. This shear may significantly impact wave dynamics, in addition to causing viscous 
heating and mixing. 



Tidal Heating in White Dwarf Binaries 23 



Puller, J., Lai, D., 2012, ApJL, 756, 17 (Paper III) 
Gilfanov, M., Bogdan, A., 2010, Nature, 463, 924 
Goldreich, P., Nicholson, P., 1989, ApJ, 342, 1079 
Golreich, P., & Soter, S. 1966, Icarus, 5, 375 
Gonzalez Hernandez, J., et al. 2012, Nature, 489, 533 

Hermes, J. J., Kilic, M., Brown, W., Winget, D., Allende Prieto, C, Gianninas, A., Mukadam, A., Cabrera-Lavers, S., Kenyon, S., 2012, 

ApJL, 757, 21 
Hut, P. 1981, A&A, 99, 126 
Iben, I., Tutukov, A. 1984, ApJS, 54, 335 
Iben, I., Tutukov, A., Fedorova, A. 1998, ApJ, 503, 344 

Kilic, M., Brown, W.R., Kenyon, S.J., Allende Prieto, C, Andrews, J., Kleinman, S.J., Winget, K.I., Winget, D.E., Hermes, J.J. 2011, 
MNRASL, 413, 101 

Kilic, M., Brown, W.R., Hermes, J.J., Allende Prieto, C, Kenyon, S.J., Winget, D.E., Winget, K.I., 2011, MNRASL, 418, L157 

Kilic, M. Brown, W., Allende Prieto, C, Kenyon, S., Heinke, C, Ageros, M., Kleinman, S., 2012, ApJ, 751, 141 

Kulkarni, S.R., van Kerkwijk, M.H., 2010, ApJ, 719, 1123 

Lorcn-Aguilar, P., Isern, J., Garcia-Berro, E. 2009, AA, 500, 1193 

Maoz, D., Sharon, K., Gal- Yam, A. 2010, ApJ, 722, 1979 

Marsh, T. 2011, CQGra, 28, 094019 

Mullally, F., Badenes, C, Thompson, S.E., Lupton, R. 2009, ApJ, 707, L51 
Ogilvie, G.I., Lin, D.N.C. 2007, ApJ, 661, 1180 

Paxton, B., Bildsten, L., Dotter, A., Herwig, F., Lasaffre, P. Timmes, F. 2011, ApJS, 192, 3 
Piro, T. 2011, ApJ, 740, L53 

Raskin, C, Scannapieco, E., Fryer, C, Rockefeller, G., Timmes, F. 2012, ApJ, 746, 62 

Schaefer, B., Pagnotta, A., 2012, Nature, 481, 164 

Segretain, L., Chabrier, G., Mochkovitch, R. 1997, ApJ, 481, 355 

Spruit, H.C. 2002, A&A, 381, 923 

Steinfadt, J., Kaplan, D.L., Shporer, A., Bildsten, L., Howell, S.B. 2010, ApJ, 716, L146 

Steinfadt, J., Bildsten, L., Arras, P. 2010, ApJ, 718, 441 

Valsecchi, F., Farr, W., Willems, B., Kalogera, V., 2012, |arXiv:1210.5023| 

Van Kerkwijk, M.H., Chang, P., Justham, S. 2010, ApJL, 722, 157 

Webbink, R.F. 1984, ApJ, 277, 355 

Willems, B., Deloye, C.J., Kalogera, V. 2010, ApJ, 713, 239 

Yoon, S.-C, Podsiadlowski, Ph., Rosswog, S., 2007, MNRAS, 380, 933 



